Microscopic Evolution of a Weakly Interacting Homogeneous Bose Gas 
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We provide a detailed description of the nonequilibrium 
time evolution of an interacting homogeneous Bose-Einstein 
condensate. We use a nonperturbative in-medium quantum 
field theory approach as a microscopic model for the Bose gas. 
The real-time dynamics of the condensate is encoded in a set 
of self-consistent equations which corresponds to an infinite 
sum of ladder Feynman diagrams. The crucial role played by 
the interaction between fluctuations for the instability gener- 
ation is thoroughly described. 
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I. INTRODUCTION 

Bose-Einstein condensation (BEG) is a remarkable 
phenomenon that has been present amongst the whole 
spectrum of research in theoretical physics since its pro- 
posal . It has challenged experimentalists with its eva- 
sive nature until about five years ago, when it was fi- 
nally observed in a series of elaborate experiments with 
weakly-interacting atomic vapours confined in magnetic 
traps and cooled down to the realm of temperatures of 
the order of fractions of microkelvins |^ . Such an aston- 
ishing environment was made possible by the combina- 
tion of refined techniques such as laser cooling, magneto- 
optical trapping and evaporative cooling. Since this re- 
markable experimental accomplishment, a great amount 
of theoretical investigation has been stimulated (see ||] 
for a comprehensive review). In part, this interest grew 
up because several experimental features can be deter- 
mined by fine-tuning various interesting parameters with 
a high level of control and accuracy. Moreover, if one 
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really intends to describe the actual evolution of the con- 
densate formation, then finite-density, non-zero temper- 
ature, and nonequilibrium dynamics effects will have to 
be taken into consideration. This may be accomplished 
quite naturally by in-medium nonequilibrium quantum 
field theory methods. This set of characteristics makes 
BEG one of the most attractive and promising systems in 
which one can use models and approximations that could 
also prove useful in very different environments such as 
neutron stars or heavy- ion collisions Q| . 

Recent experiments with dilute atomic gases were 
able to start probing quantities which are relevant to the 
understanding of the underlying dynamics of BEG, such 
as the time scales for the condensate formation and its 
final size. On the theoretical side, the microscopic be- 
havior of such an environment can be appropriately de- 
scribed by the nonequilibrium Schwinger-Keldysh closed 
time-path formalism in the quantum field theory ap- 
proach [^,0. In fact, the first steps in this direction were 
performed by Stoof and provided a qualitative idea 
of the various times scales involved in the BEG process. 
However, a microscopic approach to the condensate onset 
and growth which can explain the mechanism of instabil- 
ity generation and describe the time evolution that fol- 
lows in a quantitative way is still crude and incomplete, 
even in the case of a homogeneous gas. The second stage 
in the whole process, i.e. the growth of the condensate, 
was considered by Gardiner et al. ^ through a quantum 
kinetic formalism and the construction of a master equa- 
tion for a density operator describing the state of the 
condensate, which is equivalent to a Boltzmann equation 
describing a quasi-equilibrium growth of the condensate. 
Nevertheless, as pointed out by Stoof the simplest 
formulations based on kinetic theory do not allow for 
the observation of a macroscopic occupation of the one- 
particle ground state, and the question of the instability 
of the Bose gas system in the homogeneous case is a non- 
trivial one. 

In this article, we undertake the task of providing the 
simplest yet sensible nonequilibrium in-medium quan- 
tum field theory description of the dynamics of conden- 
sation of an interacting homogeneous Bose-Einstein gas 
that is abruptly quenched from a temperature far above 
to a temperature far below the critical one 0. This 
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approach yields a nonperturbative set of self-consistent 
equations which correspond to a resummation of the lad- 
der Feynman diagrams discussed in ref. In fact, this 
sum can be interpreted in a simpler, non-diagrammatic 
way by a careful analysis of the role played by the fluctu- 
ations, which shows that their interaction are crucial to 
the mechanism of instability generation in a model with 
no instabilities at the mean field level. Moreover, as will 
be described below, one can separate the fc-space into 
three regions that differ in their stability properties, the 
scale being determined by the interaction strength and 
the density of condensate. 

The article is organized as follows: Section II presents 
the microscopic model Lagrangian density and an effec- 
tive theory for the fluctuations. Section III contains the 
self-consistent integro-differential equations for the con- 
densate evolution that follow from the nonequilibrium 
approach. Besides, there is a detailed discussion of the 
assumptions and approximations that one has to make 
in order to provide an analytic treatment of the problem. 
In Section IV, we discuss our results for the density of 
condensate as a function of time and present our final 
comments. 



II. EFFECTIVE LAGRANGIAN DENSITY FOR 
THE FLUCTUATIONS 

The interaction Hamiltonian that describes a homoge- 
neous gas of interacting bosons of mass m has the fol- 
lowing second-quantized form (throughout this paper we 
use units such that fi = 1): 

(1) 

where (j){x,t) and (/)^(a?, t) are the boson anihilation and 
creation Heisenberg field operators, and V{x — x') is the 
two-body interatomic potential. Since we will consider a 
weakly-interacting dilute gas, we are allowed to use the 
fact that, in this case, only binary collisions at low en- 
ergy are relevant, so that we can approximate the actual 
interatomic potential by the much simpler form 

V{x-x')=g5{x-x') , (2) 

where the coupling constant g is related to the s-wave 
scattering lengh a through 



Then, for a weakly-interacting dilute and cold gas, we end 
up with a hard core interaction potential. In fact, one is 
reasonably safe to use such approximation since typical 
values for the dimensionless parameter which controls the 
validity of the dilute-gas approximation, the number of 



particles in a scattering volume, i.e. n|a|^, is always less 
than ~ 10~^. Inspired by the discussion above, one can 
write a field theory model for this system. It corresponds 
to the simplest model for a nonrelativistic complex Bose 
field, and its Lagrangian density is given by 

£ = <^*(^| + ^v^)^+,r^-f(r0f , (4) 

where the complex scalar field 0(x, t) represents charged 
spinless bosons of mass m, and g is the coupling constant 
defined above. In (Q) we have also explicitly introduced 
a chemical potential that guarantees a constant total 
density of particles 

(0*0) =n . (5) 

We also assume that the system is coupled to a heat bath 
environment with inverse temperature (3 = l/^ksT). 
The role played by temperature is very vague at this 
point. It will be discussed in detail later. For the mo- 
ment, it is here just to remind us that we will consider the 
nonequilibrium evolution of the system at finite density 
and nonzero temperature. 

Since we are interested in describing the onset and the 
time evolution of the Bose-Einstein condensate, it is con- 
venient to decompose the original fields </> and 4>* into a 
condensate (zero-momentum mode) part given by (fo and 
(Pq, and a fluctuation part, Lp and (p* , that accounts for 
the atoms outside the condensate. Assuming a homoge- 
neous condensate, we define the following decomposition: 

(f>{x,t) ^ (po{t) + (p{x,t) , 
rix,t)^ip*it)+ip*ix,t) . (6) 

Note that we make a time-dependent shift, i.e., we take 
(po{t) as an arbitrary function of time that will be deter- 
mined by the dynamics of the system. This procedure is 
different from the usual equilibrium analysis of systems 
that feature spontaneous symmetry breaking, where one 
usually performs a constant shift. Since we intend to 
describe the nonequilibrium evolution of the condensate, 
we must keep the time dependence. 

The substitution of the shifted fields (^) in the La- 
grangian (^ upto second order yields the Bogoliubov 
approximation for quasiparticles, which neglects the in- 
teraction between fluctuations and is reasonable only at 
temperatures well below the critical temperature for the 
condensate formation. At this level of approximation, 
the system is absolutely stable, and this corresponds to 
the stationary equilibrium situation that is achieved af- 
ter the whole dynamical evolution of the condensate has 
taken place. In order to allow for the appearance of an 
instability, which will be responsible for the onset and 
subsequent evolution of the condensate as a function of 
time, we must go beyond the Bogoliubov approximation. 
The simplest extension is to implement a mean field ap- 
proximation in the interactions between the fluctuation 
fields. It turns out that even this simple improvement in 
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our approach is able to trigger the instabihty that will 
drive the condensate formation. The mean field approxi- 
mation corresponds to a self-consistent Hartree approach, 
and amounts to the following modification in the inter- 
action term for fluctuations: 

g{Lp*ipY = Ag{ip*ip)ip*if -f g{ip*ip*)(pip + g{ipip)ip* tp* + 

- 9{'P*'P*)W- 9{w)f*V*] , (7) 

where the first terms in the rhs is taken as part of the 
quadratic Lagrangian for fluctuations, and the terms in- 
side the square brackets is taken as part of the interac- 
tion Lagrangian. The terms {(fif) and {ip*ip*) are usu- 
ally known as the anomalous density terms while {^p*ip) 
is the non-condensate density. In the following we work 
with the approximation of neglecting the anomalous den- 
sity terms in which is equivalent to the the Hartree- 
Fock-Bogoliubov-Popov (HFBP) (or dynamical Popov) 
approximation for dilute Bose gases, which is generi- 
cally valid for T <C Tc, when the contributions from the 
anomalous and non-condensate densities are very small 
as compared to the condensate density. This approach 
has also shown to give a reasonable first approximation 
for the description of the excitation spectrum at all tem- 
peratures provided the non-condensate density is calcu- 
lated self-consistently 11 1. 

From the decomposition of the fields and ^ in the 
HFBP approximation, the quadratic part of the La- 
grangian density for the fluctuations, Lq(}p,ip*), may 
then be written as 



dt 



1 

2m 
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(8) 



Here, we have used the fact that, under the field de- 
composition in the condensate and out of the condensate 
modes, the density constraint then becomes 



(9) 



Additionally, assuming that at the initial time the system 
is mostly composed of particles outside the condensate, 
~ {f*f>) (at t = 0), simple relations involving the 
generating functional for the correlation functions (see, 
for instance, the last section of chapter 2 in ||l^) allow 
us to write the total number density n of particles in 
terms of the chemical potential /i, valid in the mean-field 
approximation for the potential, as: = 2gn. This ex- 
pression does not depend on any equilibrium property of 
the system, depending only on the mean-field approxima- 
tion. Notice that this value for /j coincides with the usual 
one obtained in the equilibrium treatments of dilute Bose 
gases in the HFBP approximation 111 
isfies the Hugenholtz-Pines relation 
obtained in the equilibrium problem. In spite of this, our 
derivation is valid for the out-of-equilibrium regime and 
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should not to be confused with the equilibrium treat- 
ments. These considerations lead to the quadratic La- 
grangian for the fluctuations shown above (Eq. (0)). 



III. THE QUENCH SCENARIO, GREEN 
FUNCTIONS AND SELF-CONSISTENT 
EQUATIONS 

The scenario proposed here assumes that, for t < 0, the 
initial state is in equilibrium at a temperature Ti ^ Tc- 
At t = the system is then abruptly quenched to a much 
lower temperature Tf <^ Tc, where Tf is the temperature 
of the thermal bath in which the system is immersed (the 
one associated with the quantity (3 that appears in our 
description) and, of course, it will be the equilibrium tem- 
perature which the system will reach asymptotically as 
t 00. We should remark, then, that the initial tem- 
perature itself is not important in our study and it will 
not appear in our equations. It is only a conceptual tool 
that will help us in the definition of the initial conditions 
for the number densities, as will be clear below. In fact, 
since our approach is strongly out of equilibrium, it is not 
possible to define a "temperature parameter" when the 
condensate is evolving as it is usually done in a quasi- 
equilibrium treatment. Technically, this means that the 
final equilibrium temperature, Tf — l/{kBP), will show 
up in the nonequilibrium propagators through boundary 
conditions (see below). 

This kind of quench is easily attained in the exper- 
iments of Bose-Einstein condensation of atomic gases, 
where the typical relaxation time scales are long enough 
(around ~ 0.1s, depending on the temperature [0) to 
allow for a fast drop in the temperature of the system that 
evolves afterwards out of equilibrium. With this choice 
of initial state, it is reasonable to consider the dynamics 
of the building up of the condensate state, which at the 
initial time is ncond.(i = 0) = \'Po{t = 0)p ~ 0, and the 
depletion of the excited states, which at t = is given 
by ricxcX^ = 0) = {if* if) ~ n, as essentially a two-level 
problem. It is clear that this approximation breaks down 
for temperatures close to the critical temperature, where 
the detailed treatment would require a thorough study 
of the dynamics among the many levels of excited states. 
Furthermore, at this point the mean field approximation 
would not be reliable any more. 

In the scenario described above, the condensate builds 
up subject to the density constraint relation, which may 
be expressed in terms of the averages of the real and 
the imaginary parts of if and tp* ((p = ipi + ip2 and 
p* = Pi — 1^2)- Spatial translation invariance yields: 



Xt) 



noxc.(0 = (pi{t)pi{t))^„ + {p2{t)ip2{t))ipo 



(10) 



In order to express these self-consistent equations in 
a concrete and convenient form let us define the field 
averages in terms of the Green functions for the fields. 
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The Green functions can be defined from a generating 
functional for our Lagrangian model, Eq. (jj) , where the 
generating functional Z[J], in terms of an external source 
J, is given by 



Z[J] 



D(f>D(f>*e^p{iS[(l),(l)*,J]} , 



(11) 



where the classical action is given by 



^[(/),r,J]= / d^x{C[(l),(j)*] + Jix)(b*ix) + J*{x)(b{x)}. 
Jc 

(12) 

In the time integration is along a contour suitable 
for real-time evaluations, which we choose as being the 
Schwinger closed-time path ^,0], of common use in con- 
densed matter non-equilibrium problems, where the time 
path C goes from — oo to -l-oo and then back to — oo. The 
functional integration in (^l]) is over fields along this time 
contour. As in the Euclidean time formulation, the scalar 
field is still periodic in time, but now (j){t, x) = <l)[t—il3, x). 
Temperature appears due to the boundary conditions, 
but now real time is explicitly present in the integra- 
tion contour. Denoting by J+ and J" the sources in the 
— cx) to +00 path and -\-oo to —00 paths, respectively, 
along the time path contour, the generating functional 
Z[J+, J~] can be written as 



D(t)D(t)* exp^i j d'^x{C[(j)+,(j}X,J+] 
01, J-]} ) . (13) 



Z[J+,J- 



By performing the path integrals over the quadratic 
forms as usual, one obtains that the generating func- 
tional in the Schwinger's closed-time path formalism can 
be written as (with 0, (jf written in terms of the real 
fields 01 and 02) 



Z[J+,J"]=cxpi i j d'^ 



L 



int I ^ 




SJ- 



X exp - / d^x d\jnx)GfAx, y)j'Ay) (14) 



the periodicity condition in real time, Kubo-Martin- 
Schwinger (KMS) condition. Here, /3 is the inverse of 
the temperature of the thermal bath and appears as a 
consequence of the boundary conditions arising from the 
construction of the complex time path. 

The field averages of last section can then be expressed 
in terms of the Green functions for ipi and (^2 as (j = 1, 2) 



{ipj{t)LPj{t)) 



/(03HG>(k,M) 



(16) 



and are constructed from the homogeneous 
solutions to the operator of quadratic fluctuation modes 
which, using Eq. (^ expressed in terms of tpi and ip2, 
are given by (in momentum space) 



dxi(k,t) 



dt 



— -g\^o\^]x2(.Kt) = 



(17) 



The boundary conditions for the solutions of the equa- 
tions above are such that, for t < 0, |(/3o(OP = 
and xi(k, i) and X2(k, i) can, for example, be given by 



Xi(k,t) =e-x7V2iVand X2(k,i) 



le 



/y'2N, where 



Ek = k^/(2r7i) and the normalization factor N can be 
fixed by imposing the constraint condition, Eq. (|^), at 
t = ( |((5o(0)p = 0). We should stress that, from our 
choice of initial configuration, we have a well-defined ini- 
tial state, given by no ~ and Uexc — '^■j and a well- 
defined equilibrium final state, at temperature Tf -^Tc, 
towards which the system will evolve. The final situation, 
then, is an equilibrium one, which allows for the defini- 
tion of a final equilibrium density matrix, p — Tr e~^-^ , 
where (3 is the inverse of the temperature of the heat 
bath that was responsible for the quench {Tf <C Tc). 
Therefore, one can consistently define the Schwinger- 
Keldysh closed time-path procedure in order to describe 
the noncquilibrium evolution that will follow. The condi- 
tions fulfilled by the noncquilibrium propagators in that 
formalism are analogous to the KMS conditions obtained 
in the imaginary time description of equilibrium prob- 
lems. 

By decoupling the set of equations in (p7[), we obtain 



with a,b = +, — and j = 1, 2. The Green functions that 
enter the integrals along the closed-time path contour in 
( |r^ ) are given by (in momentum space) 

G++(k,t,i') = G>{k,t,t')0{t - t') + G< (k,i,t')0(t' - t) 
G;r (k, t,t')^ G> (k, t, t')e{t' - i) + G< (k, t,t')9{t-t') 



G+-(k,t,t')--G<(k,i,t'), 



33 ^ 



G-/(k,t,i') = -G>(k,i,<') 



(15) 



The functions G^ and G^ satisfy the property 



G<(k,i,i') 



G^ (k, t — i(3,t'), which is recognized as 



Xi 

X2 



— Xi +^+uJ-Xi = 



X2 +UJ+UJ-X2 







(18) 



where uj± — k^/(2m) ± g|yop. Defining ^1 = xi/y/^ 
and ^2 = X2/yA^, Eq. (18) can also be written as 



6 + Uj1^2 



(19) 



where 
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-'1(2) 



2uj_ 



(+) 



(+) 



(+) 



(20) 



We can see from these decoupled equations that those 
modes with k^/(2m) < gl^/JoP are the unstable ones and 
they will drive initially the excited particles towards con- 
densation. At this point, one should recall that the in- 
stability arises as a consequence of the interaction be- 
tween fluctuations. Should one discard those interac- 
tions, the mean- field description would result in a abso- 
lutely stable system. This situation is completely differ- 
ent from the one encountered in the relativistic approach 
| p5[ . Those unstable modes, for which k"^ / {2m) is smaller 
than (7|(/7op, correspond to the (exponential) growth of 
the long-wavelength fluctuations which drive the process 
of phase transition, or condensation in our case. This is 
similar to the phenomenon of spinodal decomposition in 
statistical mechanics, typical of second order phase tran- 
sitions [ pj[ |. 

In terms of the mode functions xi and X2^ and tak- 
ing into account the boundary conditions for them, the 
Green's functions can be expressed as (j = 1, 2) 



G>(k,i,t') 



1 — 

e-P^-X*{^.t)xj{^X)\ 



(21) 



The Green's functions for the complex fields ip* are 
expressed as usual by 

^[-zG>^.(k,i,t)] , (22) 



{V*{t)^{t)) 



cPk 
(2^ 



-iG> ^(k,t,i)] . (23) 



In terms of (|2^) and (|2^), we have that the non- 
condensate density noxc(i) can be written as 



"-oxc(0 = {^Pl{t)Lpi{t)) + ((p2{t)ip2it)) = 

= i((^(%*(t)) + (^*(%(t))) 



(24) 



and the KMS condition can be expressed, in this case, 
as G>^,ik,t - i(3,t') = [G> ^(k,t,t')]*, or G>.^(k,t - 

j/3,i') = [G>^.(k,i,t')]*- 

Using Eq. (pl|), we then obtain that 7ioxc(^) can be 
expressed as 



^exc (0 



\xiiKt)\^+ 



(27r)3 

X2{Kt)f] coth(^) 



(25) 



At t 0, from the initial conditions imposed over the 
system, the system is all in the non-condensate state and 



the condensate density Uc 



ml 



0; xi and X2 are 



nexc(O) ^ n = (v^V) ^ (V'?) 1,^0 + i^Dl^, , (26) 
which, from Eq. (p5|), gives 
1 



27V7r2 



dk fc2coth(^^^ 



(27) 



The T = part of (^) is divergent and represents the 
zero-point energy. By subtracting the zero-point energy 
we obtain for the normalization N 



(28) 



Using the expression for the density in terms of the 
critical temperature for condensation of an ideal Bose 

gas (l//3o) 0, n = C(3/2) [m/{27TPo)f\ we find a 
simple expression for the normalization N as given by 
N ~ 2{Po/P)^, which corresponds to the pre- factor 
used in Ref. Using these constraint conditions to- 

gether with the initial boundary conditions, from Eq. 
(Eq), we can then express the condensate density as 



1 f ry/ST^alvom 



purely oscillatory functions. We then obtain that 



X [l-iV(|Xi(k,t)|2-K|x2(k,t)|2)]nk(/3) , (29) 

where nk(/3) — [e^'^^ — 1)^^- The ultraviolet divergences 
were explicitly cancelled by subtracting the T = com- 
ponent of (|25|). This is fine, since the contributions from 
the non-condensate part of the density are neglegible for 
weak interacting dilute Bose gas systems at T = [1). 
Note also that, in our out-of-equilibrium approach, there 
are no infrared divergences since the finite time is a nat- 
ural regulator. However, for the equilibrium t — *■ c», the 
critical temperature will be modified by the interactions 
as pointed out in Ref. ||l9|,^ . 

Equations (|l^) and (^9|) form an integro-differential 
system that may be solved for (^o(0 numerically, given 
the initial conditions for ipo{t), XiO^i and X2(k, t) men- 
tioned before. Indeed, this system of equations deter- 
mines completely the time evolution of the condensate 
density as a function of the temperature and of the total 
density of the gas. At this point, one could ask about 
the contribution from those terms responsible for the de- 
cay of quasiparticles which are related, microscopically, 
to the imaginary part of the self-energies in higher or- 
der perturbation theory. These contributions are impor- 
tant to describe the detailed dynamics for temperatures 
close to the critical temperature, which is not the case 
we have at hand. In such range of temperatures, the 
typical time scales for condensation are long enough so 
that scattering and growth among adjacent levels, which 
are described by those microscopic processes, are of rel- 
evance to the dynamics in this regime. Those processes 
have been described in the literature through the use of 
quantum kinetic theory, i.e. quantum Boltzmann equa- 
tions, and could as well be done within the approach 
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described here when considering higher-order corrections 
to our field averages. 

Equation ( p5| ) is the first order term in the finite- 
temperature quantum field perturbation expansion for 
((/?|). Higher-order corrections for the equal-time two- 
point field averages can be expressed in terms of the co- 
incidence limit of the (causal) two-point Green functions 
G^pip' and G^p'^, which satisfy the Dyson equation (the 
indices stand for ip and ip*): 



(30) 



where is the (matrix) self-energy, and G^^ is the 
zeroth-order non-interacting Green function, satisfying 
the equation (in momentum space) 



±i- 



dt 



{k,t,t') = S{t~t) 



(31) 



One of the advantages of expressing the Green's func- 
tions in terms of the solutions of (jT^) is the possibility 
of obtaining, in an unambiguous way, all higher-order 
corrections to the two-point and many-point functions. 



IV. RESULTS AND DISCUSSION 

Prom Eq. (|2^) we find that initially the dynamics 
will be dominated by essentially two regimes in the k- 
space: for {k'^/2m) <^ ffl^'oP, or spinodal regime, is the 
regime caracterized by the exponential growth of the fluc- 
tuations, driving the system towards Bose- Einstein con- 
densation. Following a crossover at [k"^ /2m) ~ 5|<PoPj 
the modes for which (fc^/2m) ^ 5|<PoP will be oscilla- 
tory stable modes. The energy of these modes are high 
enough such that they effectively sample a symmetric po- 
tential and will not contribute to the dynamics. These 
properties determined our choice of momentum cut-off 
in Eq. (p9|). Equations ( p7| ) and (^) form a set of 
integro-differential equations that can be solved numer- 
ically, given the initial conditions mentioned before. In 
Fig. 1 we show the results for the condensate density 
|</3oP, as a function of time, for different values of tem- 
perature. From these results we notice the characteris- 
tic growth curves for the condensate. For very initial 
times, the system is dominated by the stable modes. As 
the condensate develops the number of unstables modes 
increases. At some point they dominate and make the 
condensate grow exponentially fast, until the equilibrium 
state is reached. This state is determined by the temper- 
ature at which the system is quenched below the critical 
temperature, and by the total density. 

It is important to point out that the evolution of the 
condensate is completely driven by the interaction be- 
tween the microscopic fluctuations of the field around 
the condensate. In fact, in this non-relativistic system, 
the interaction between fluctuations is responsible for the 
appearance of the instability which allows for the onset 



of the condensate (see, e.g., for a description of what 
happens in the relativistic case). One should also note 
that not all the excited particles condense, since there 
will always be a fraction (which depends on various pa- 
rameters for a particular system, and on the temperature 
of the thermal bath) of excited modes, with high enough 
frequency, that remains stable. 
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FIG. 1. Condensate density as a function of time for 
na^ = 0.01 and Ti/Tq = 0.06, T2/T0 = 0.08 and T3/T0 = 0.1. 
Here, t = {h/ma?)t is a dimensionless time and po = ci''|(^oP 
is a dimensionless density. 

The system of equations we have derived is completely 
nonperturbative and able to describe the time evolution 
of BEG in the extremelly far from equilibrium regime, 
where standard methods based on the mean-field Gross- 
Pitaevski equations are not applicable. In this sense, our 
description should complement the usual approach via 
Boltzmann equation. Furthermore, the system of equa- 
tions obtained here, can be viewed as a resummation of 
the ladder Feynman diagrams mentioned by Stoof in . 

In summary, we have developed an out-of-equilibrium 
non-perturbative quantum field theory description of 
the condensation process of an interacting homogeneous 
Bose-Einstein gas quenched below the critical tempera- 
ture. Although we focused this article on the instabil- 
ity process that generates the condensate (i.e., the short 
time behavior), for t — > cx) our results confirm the ex- 
pected behavior for this limit However, the equilib- 
rium (i — > 00) values of the condensate fraction are lower 
than the experimental results jl^ and the calculations of 
Dalfovo et. al. [^j. This may be due to our approxima- 
tion of neglecting incoherent coUisional processes, which 
is a valid approximation in an infinite homogeneous gas 
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at very low temperatures and densities, but otherwise 
may give an important contribution. We expect that the 
self-consistent inclusion of pair terms should account for 
most of these contributions. 

In spite of the absence of non-homogeneity effects, we 
hope that the approach developed here may be useful in 
the analysis of transients in realistic Bose-Einstein con- 
densation experiments with atomic gases. Our results 
should also be applicable to trapped atomic gases in the 
central region of wide traps. Moreover, with a suitable 
generalization of the formalism presented above, we could 
be able to develop a theoretical description of the dy- 
namical aspects of a recently proposed experiment [2l|, 
regarding the Bose-Einstein condensation in a weakly- 
interacting photon gas in a nonlinear Fabry-Perot cavity 
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